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We consider the interaction of a quantum mechanical SQUID ring with a classical resonator (a 
parallel LC tank circuit). In our model we assume that the evolution of the ring maintains its 
quantum mechanical nature, even though the circuit to which it is coupled is treated classically. 
We show that when the SQUID ring is driven by a classical monochromatic microwave source, 
energy can be transferred between this input and the tank circuit, even when the frequency ratio 
between them is very large. Essentially, these calculations deal with the coupling between a single 
macroscopic quantum object (the SQUID ring) and a classical circuit measurement device where 
due account is taken of the non-perturbative behaviour of the ring and the concomitant non-linear 
interaction of the ring with this device. 



INTRODUCTION 



With the now very serious interest being taken in the 
possibilities of creating quantum technologies such as 
quantum information processing and quantum comput- 
ing [1-3], much attention is being focused on the applica- 
tion of Josephson effect devices, particularly the SQUID 
ring. As has been demonstrated recently, with the appro- 
priate ring circuit parameters and operating temperature 
these devices can display manifestly macroscopic quan- 
tum behaviour such as macroscopically distinct superpo- 
sition of states [4-11]. In any considerations of quan- 
tum technologies the role of the environment, as coupled 
to the quantum object, has featured very strongly [12- 
14]. With regard to using superconducting systems in 
quantum technologies, it has been shown that Josephson 
weak link circuits, and in particular SQUID rings in the 
quantum regime, are highly non-perturbative in nature 
and can generate very strong non-linear interactions with 
classical circuit environments [15]. In this paper we pro- 
vided a demonstration that this non-perturbative (non- 
linear) behaviour is crucial to the understanding of the 
interaction of SQUID rings with circuit environments. In 
this work we first consider the adiabatic (ground state) 
interaction of a quantum regime SQUID ring inductively 
coupled to a classical parallel resonance LC (tank) cir- 
cuit. This circuit system is shown schematically in fig- 
ure 1(a). We then continue to consider energy trans- 
fer between a classical monochromatic microwave source 
(the input) and this classical electromagnetic (em) field 
mode (the output) via the quantum mechanical SQUID 
ring (figure 1(b)). Here, as in other work [16], we have 
modelled this output mode as a parallel LC equivalent 
circuit. However, in this paper we have treated this cir- 
cuit environment as classical since we wish to explore 
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FIG. 1: (a) Schematic of a (quantum regime) SQUID ring in- 
ductively coupled to a driven tank circuit (two mode system) 
(b) schematic of a three mode system: classical input mi- 
crowave excitation source - SQUID ring - and an inductively 
coupled tank circuit. 



a region of the system parameter space, involving large 
input-output frequency ratios, which is of current exper- 
imental significance. When treated fully quantum me- 
chanically, such simulations are beyond the reach of the 
computing power available to us. As the essential re- 
sult of this paper, we show that the non-perturbative 
(non-linear) interaction generated by the SQUID ring is 
made manifest through energy conversion between the in- 
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put and output, even if these differ greatly in frequency. 
In part this has been inspired by past experimental re- 
sults where we have demonstrated [17] that very high ra- 
tio frequency down conversion can occur between a mi- 
crowave source and a radio frequency (w 20MHz) tank 
circuit via the intermediary of a very small capacitance 
SQUID ring [17]. In these previous results we recorded 
a frequency down conversion ratio of 500:1. In recent, as 
yet to be published data, this has been improved upon 
by refinements in electronic technique to yield a ratio as 
high as 6,500:1 between a a few GHz microwave input 
field and a sub-MHz tank circuit. We know of no clas- 
sical non-linear circuit device that could generate such 
large ratios [18] and, given the very small capacitance 
(niobium point contact) SQUID rings we were using in 
our experiments, it was clearly of interest to see whether 
a quantum regime SQUID ring could lead to this mas- 
sive frequency down conversion. In a recent publica- 
tion [19] we considered this problem in the context of 
the persistent current qubit circuit model introduced by 
Orlando et al [20], where this quantum element acted as 
the medium for coupling two classical field environments 
together. Specifically, the qubit was driven by a classi- 
cal electromagnetic (em) field (the input) which was used 
to pump the qubit into an excited state. To deal with 
energy conversion to an dissipative output em mode, a 
quantum jumps approach was adopted to model the de- 
cay of the excited qubit with the energy dumped into this 
mode. Using this approach we were able to demonstrate 
small ratio (500MHz to 300MHz) frequency down con- 
version from input to output through the intermediary 
of the quantum qubit loop. In this work, with a single 
weak link SQUID ring as the quantum intermediary, we 
adopt a less complicated model, where the SQUID ring 
simply follows Schrodinger evolution without the intro- 
duction of quantum jumps. We then show, to the limits 
of the computational power available to us, that a two 
orders of magnitude frequency down conversion (from mi- 
crowave to radio frequencies) is possible. In doing so we 
also demonstrate that the interaction between a quantum 
mechanical SQUID ring and its classical circuit environ- 
ment is in no way trivial and requires a detailed under- 
standing of the non-perturbative properties of SQUID 
rings. We emphasise that with this result established, 
and without any other obvious physical constraints, it 
seems reasonable to assume from a theoretical viewpoint 
that even higher ratio frequency down conversion could 
be realised if the necessary computational power were 
available. At present, apart from essentially fixed fre- 
quency masers, sources of quantum mechanical em fields 
at microwave frequencies are not readily available. Thus, 
given the non-perturbative properties of SQUID rings in 
the quantum regime, it also seems reasonable to consider 
how such devices could make manifest their quantum me- 
chanical nature through the interaction with classical em 
fields and field oscillator modes. 
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FIG. 2: Plot of the eigenenergies of the SQUID ring Hamil- 
tonian, found by solving the TISE, as a function of external 
applied flux with hu a = 0.006$o/A s and hv = 0.035$o/A s 



In the well known lumped component model of a quan- 
tum mechanical SQUID ring (in this work one Josephson 
weak link device, of effective capacitance C s , enclosed 
by a thick superconducting ring of inductance A s , with 
a classical magnetic flux of Q x applied to the ring) the 
Hamiltonian for the ring is given by [21, 22] 



H, 



2C S 



2A S 



hv cos 



2tt$ s 



(1) 



where, with a circumflex denoting operators, Q s (the 
magnetic flux threading the ring of the SQUID device) 
and Q s (the Maxwell electric displacement flux between 
the electrodes of the weak link) are canonically conjugate 
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quantum variables such that 

—ihd/d& s and A& S AQ S > h/2. Here, the third term on 
the right hand side of (1) is due to the Josephson phase 
coherent coupling energy, with matrix clement hv/2, aris- 
ing from the weak link critical current I c = 2ev, with pe- 
riodicity set by the superconducting flux quantum $o = 
h/2e. We assume that the ambient temperature (T) of 
the SQUID is such that Hoj s 3> ksT for a character- 
istic ring oscillator frequency of lo s /2tt (= l/2n^/A s C s ) . 
Operating in this quantum regime the time independent 
Schrodinger equation (TISE) for the SQUID ring then 
reads 



2C 



2A 



2tt$ s 



\tp K ) = E K \ip K ) 
(2) 

for ring eigenfunctions \ip K )(n = the ground state, n = I 
the first excited state, etc.) and ring eigenenergies E K (as 
shown in figure 2), these eigenenergies being <I>o-pcriodic 
in the external magnetic flux applied to the ring. 
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DYNAMICS OF A COUPLED QUANTUM SQUID 
RING-CLASSICAL RESONATOR SYSTEM 



Adiabatic Regime 

In the time independent (adiabatic) case the experi- 
mentally accessible state is the ground state (k = 0) for 
which the expectation value of the screening supercurrcnt 

flowing around the ring is (I s (3? x )) 



lO&s — $ x ^j /A s \ [17, 21]. For relatively large values 

of hv (and hence I c ) this ground state screening cur- 
rent takes the form of a rounded sawtooth centred on 
&x = $0/2 (modulo n<f> , n integer). Thus, the ring 
screening current is clearly a highly non-linear function of 
$x, as is its magnetic susceptibility Xs = d(I s (& x )) /d& x . 
At <£>£ = (n + 1/2) $0, n integer, the ring is in an quan- 
tum superposition, with equal coefficients, of clockwise 
and anti-clockwise screening currents (equivalently flux 
states). A rounded sawtooth jump in the screening cur- 
rent, centred at half integer bias flux, correspondingly 
generates a narrow positive spike in Xs (®x) at this flux 
value. At integer bias flux (& x = n& ) the susceptibil- 
ity is a minimum, increasing monotonically to a maxi- 
mum at $ T = (n+ 1/2) $o- As we have shown, in this 
simple ground state model of a quantum regime SQUID 
ring the evolution of its quantum state can be monitored 
through its effect on the classical dynamics of a mea- 
surement circuit coupled to it [15, 17, 21]. Typically 
this takes the form of a low frequency, parallel LC tank 
circuit inductively coupled to the SQUID ring and ex- 
cited at, or very close to, its resonant frequency. The 
resonant frequency of this tank circuit (wj/27r) is usu- 
ally extremely low compared with the natural oscillator 
frequency of the SQUID ring (with the implicit assump- 
tion that hwt <C ksT). This configuration is very well 
known and forms the basis for one type of ultra sensitive 
magnetic flux detector - the ac or radio frequency bi- 
ased SQUID magnetometer [22]. With this tank circuit 
playing the role of the classical measurement system for 
the SQUID ring, a changing Xs (&x) leads, in the ground 
state of the ring, to downward shifts (due to the positive 
Xs (®x)) of the resonant frequency of the tank circuit. 
For this ground state, therefore, following the change in 
the frequency (and amplitude) of the tank circuit reso- 
nance allows us to extract some information concerning 
the evolution of the quantum state of the ring (i.e. the 
coefficients in its screening current superposition) as a 
function of <f> T . 



(dH/d$ x 



often convenient to introduce an approach well known 
in atomic physics. These capacitances can be consid- 
ered as measures of the effective mass of each circuit 
(SQUID ring and tank circuit resonator). Crucially, for 
quantum regime SQUID rings these capacitances (effec- 
tive masses) will be markedly different in magnitude, e.g. 
for the ring typically a few xlO~ 15 F and for the tank 
circuit sa a few xlCU n F, or more. This low mass-high 
mass (here, low capacitance-high capacitance) situation 
was dealt with many years ago by Born and Oppenhcimcr 
for nuclear-electron motion. Adopting this approach we 
compute from solutions of (2) the expectation value of the 
screening supercurrent in the SQUID ring (fast) and sub- 
stitute this as a feedback term into the classical equation 
of motion for the tank circuit (slow). Provided that the 
resonant frequency of the tank circuit is low enough, so 
that to an extremely good approximation the quantum 
SQUID ring remains adiabatically in its ground state, 
the equation of motion for the coupled ring-tank circuit 
system then reads [15, 21] 



Ct 



d 2 $ t 
dt 2 



1 d<& t 



1 



Rt dt L t 



+ —<5>t=hn{t)+llU Is V) (3) 



where L t and Ct are, respectively, the tank circuit in- 
ductance and capacitance, $4 is the magnetic flux in the 
tank circuit inductor and Rt is the resistance of the tank 
circuit on resonance. In this Born-Oppcnhcimcr approx- 
imation we assume that the ring remains in one of its 
instantaneous energy eigenstates, i.e. solutions of (2); in 
reality this means the lowest energy state (the adiabatic 
limit). As we have demonstrated in the literature [23], 
this approximation holds very well if the frequency of the 
drive is very small indeed compared with any separations, 
in frequency, between the SQUID ring energy levels. This 
certainly appears to be the case for radio frequency (rf) 
circuits, typically used to probe single weak link SQUID 
rings, as well as in their application in SQUID magne- 
tometry [24, 25]. We note that since the SQUID ring is 
macroscopic in nature (as, of course, is the tank circuit) 
there exists a significant back reaction between these two 

circuits, as evidenced by the fi lip I s term in (3). In 

practise this means that the SQUID ring (through the co- 
sine in its Hamiltonian) generates a non-linear dynamic 
in the classical circuit environment to which it is coupled, 
in this case a tank circuit. For small tank circuit drive 
amplitudes this manifests itself as a frequency shift in the 
power spectrum of the tank circuit as shown in figure 3 



The Born-Oppenheimer Approximation 

Given that the oscillator frequencies of both the 
SQUID ring and the (classical) tank circuit resonator are 
functions of their circuit capacitances, it is, in this limit, 



TIME DEPENDENT SCHRODINGER 
EQUATION DESCRIPTION - EXCHANGE 
REGIONS 

Following on from the time independent Schrodinger 
description of a SQUID ring (2), the time dependent 
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FIG. 3: Frequency shift of a 130MHz tank circuit (Lt c — 
3 x 10~ a H, Ctc = 5 x KT^F), coupled to the SQUID ring 
of figure 2, as a function of external applied flux for the ring 
in its ground state only (figure 2(a)) calculated using the the 
Born-Oppenheimer (adiabatic) approach. 



FIG. 4: Time averaged energy expectation values of the ring 
in the presence of a microwave source (input) at a frequency 
of 13GHz using the first two energy eigenstates as initial con- 
ditions. 



Schrodinger equation (TDSE) for the ring takes the form 




where <fr x = $^, tat + $ (t) is now comprises of a time de- 
pendent magnetic flux $ (t) together with a static compo- 
nent $!j. tat . The intrinsically extremely non-pcrturbativc 
quantum nature of the SQUID ring would indicate that 
the application of cm fields low in frequency compared 
with the frequency difference (xh) between adjacent en- 
ergy levels in the ring should still induce transitions be- 
tween its quantum levels [26]. Thus for level differences 
of a few hundred GHz (quite typical of quantum regime 
SQUID rings), we expect em fields at microwave fre- 
quencies (a few GHz) to induce transitions, provided 
the field amplitude is high enough [26]. As our exam- 
ple in this work we consider a SQUID ring with circuit 
parameters commensurate with operation in the quan- 
tum rcgimc[2, 27]. For this we choose C s = 5 x 10~ 15 F, 
A s = 3 x 10~ 10 H and hv = 0.035$g/A s , the latter yield- 
ing a weak link critical current I c (= lev) of 2^iA. Fur- 
thermore, for an oxide insulator, tunnel junction weak 
link (oxide thickness w lnm, oxide dielectric constant 
« 10) in the SQUID ring, this capacitance corresponds 
to junction dimensions 0.25/im square, well within the 
capabilities of modern microfabrication technology. In 
turn, these dimensions imply a maximum supercurrent 
density in the weak link of around 4kAcm -2 , a perfectly 
reasonable value for current experiments on quantum 
regime SQUID rings. The SQUID ring circuit param- 
eters C s = 5 x 10" 15 F and A s = 3x l(r 10 H lead to a 
ring oscillator frequency lo s /2tt = 130GHz. If we assume 



that the (planar) SQUID ring is fabricated in niobium, 
which is often the case, then this ring oscillator frequency 
is small compared with the superconducting energy gap 
in this material (w 1000GHz) at low reduced tempera- 
tures T/T c , where T c is the superconducting critical tem- 
perature (9.2K for niobium) for an operating (ambient) 
temperature in experiment typically around 40mK. 

Using the ring circuit parameters given above we 
showed in figure 2 the lowest few energy eigenvalues of 
the SQUID ring Hamiltonian, found by solving the TISE 
(4), as a function of the static applied magnetic flux <&| tat . 
Although such solutions of the TISE may be adequate for 
modelling this system if the em fields frequencies are ex- 
tremely small compared with frequencies (xh) separating 
the ring energy levels (i.e. in the Born-Oppcnhcimcr ap- 
proximation, above), in general we must solve the TDSE 
where time dependent fields are involved. This is clearly 
demonstrated in figure 4. Here, using the the ring pa- 
rameters of figure 2, we display the time averaged energy 
expectation values of the TDSE (4) as a function of <J)| tat , 
for an applied microwave field of frequency 13GHz, am- 
plitude 0.001<I>o and using the first two eigenstates of the 
Hamiltonian as initial conditions. As can be seen, for 
this microwave frequency and amplitude there is signifi- 
cant energy exchange between these time averaged ener- 
gies which takes place over a a very narrow range of the 
bias flux & s J at . We term these exchange regions [16, 26], 
or equivalcntly, transition regions. We note that at the 
centres of these exchange regions the corresponding fre- 
quency differences between the energy eigenenergies in 
figure 2 are close to integer multiples of the applied mi- 
crowave frequency. At sufficiently high em frequencies 
and amplitudes several, or many, exchange regions are 
generated with an energy spacing between adjacent re- 
gions very close to hw mw , where ui mw is the (angular) fre- 
quency of the em field. As we have shown in a previous 
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publication [16], it is in these exchange regions that the 
interaction between a SQUID ring and one or more em 
fields, is significant. For example, when a field mode, or 
modes, is treated quantum mechanically, entanglement is 
generated between the ring and the mode(s), reaching (as 
with the strength of the interaction) a maximum in the 
centre of an exchange region. Given the result shown in 
figure 4, where it is clear that energy is being exchanged 
between the SQUID ring and the (classical) em field, we 
now consider whether a similar exchange can take place 
between this input field and a classical output field mode, 
coupled together through the ring. 



this experimental background [17], we consider the in- 
teraction and energy transfer via the SQUID ring be- 
tween a microwave input, acting purely as a source, and 
a lower frequency, resonant circuit output. This mod- 
els the SQUID ring-tank circuit (ac-biased) magnetome- 
ter configuration [24] in a non-adiabatic em (microwave) 
field. In this paper, we assume that there is a macroscop- 
ically significant back reaction between the tank circuit 
and the SQUID ring [15, 28]. However, for simplicity, 
we assume that there is no direct coupling between the 
cm input and output fields and note that this is easy to 
realise experimentally 



Non-adiabatic Regime 

We now investigate a scenario, accessible at the cur- 
rent level of experimental practise [17], of a classical 
field input and a classical field mode output, differing 
widely in frequency. We demonstrate that energy can be 
exchanged via a highly non-pcrturbative, quantum me- 
chanical SQUID ring. In particular, and in line with 



In order to deal in a general manner with this coupled 
ring-tank circuit system we can no longer assume that 
the behaviour of this system is adiabatic in nature. By 
implication, this means that the tank circuit resonant fre- 
quency is sufficiently large compared with the frequency 
separations between the SQUID ring levels. Hence, we 
can no longer invoke the Born-Oppcnhcimcr approxima- 
tion. The TDSE for the system then takes the form [24] 



01 

2C 



($ s - (^ 



2A fl 



— fw cos 



2tt$ s 



d\ 4>(t)) 
dt 



(5) 



where, as can be seen, the total external magnetic flux 
applied to the SQUID ring consists of static bias, mi- 
crowave excitation and back reaction tank circuit contri- 
butions, i.e. $ x = <&! tat + $™ w + ^$ t . Again, as with the 
Born-Oppcnhcimcr approximation (3), we assume that 
this back reaction is macroscopically significant. 

Given that the SQUID ring is now allowed to follow a 
Schrodingcr evolution and retain its time dependence (3) 
in this non-adiabatic regime, the equation of motion for 
the tank circuit now becomes 



d 2 $f 1 d$ f 1 



L (t) V (t) 
(6) 



where 



I. = - 



dH s = $ s - $ x jt) 
d<I> x A s 



We solve the simultaneous coupled differential equations 
(5) and (6), where the frequency ratio between input mi- 
crowave and output tank circuit modes is large (> 100). 
In our opinion these results shed light on both the gen- 
eral problem of the description of the quantum-classical 
interface and, in particular, the interaction of non-linear 



devices such as SQUID rings with their classical environ- 
ments. 



RESULTS 

As a first check of the validity of this model we com- 
pare it with the established Born-Oppcnhcimcr approx- 
imation for the limiting regime where there are no mi- 
crowaves applied and the tank circuit is driven at a fre- 
quency ft so low that, to an extremely good approxima- 
tion, the SQUID ring remains adiabatically in its ground 
state. This holds even though, through solutions of (2), 
a whole spectrum of ring cigenstates is available. In or- 
der to demonstrate the correspondence between these two 
models in the limit of low tank circuit drive frequency, we 
calculate the 4>| tat -dependent frequency shifts f (Q s x tat ) 
using the ring parameters of figure 2 and a tank circuit 
resonant frequency 130MHz and a /i = 0.01. In figure 5 
(as we did in figure 3) we show the computed ring-tank 
circuit resonant frequency shift in this limit as a func- 
tion of <f>| tat for both the Born-Oppenheimer approach 
(in blue) and our generalised model (in red). It is quite 
apparent that there is a very high degree of agreement 
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FIG. 5: Frequency shift of a 130MHz tank circuit (Lt c — 
3 x W~ S H, C tc = 5 x 10 _11 F), coupled to the SQUID ring 
of figure 2, as a function of external applied flux for the ring 
in its ground state only (figure 2) calculated using the both 
the Born-Oppenheimer (in blue) and non-adiabatic (in red) 
approaches. 



between the two approaches in this limit. As our example 
of the frequency conversion process we consider a SQUID 
ring with circuit parameters used in the computed results 
of figures 2 and 4. For the SQUID ring parameters we 
have chosen the characteristic ring oscillator frequency 
to be close to 130GHz. Referring to the computed re- 
sults of figure 4, we again choose the frequency for the 
microwave source to be 13GHz, i.e. an order of magni- 
tude lower than the ring characteristic frequency, and set 
the tank circuit resonant frequency another two orders of 
magnitude lower (i.e. 130MHz). Within the limits of the 
computational power available to us this provides us with 
the opportunity to demonstrate extreme (1000:1), quan- 
tum SQUID ring mediated, energy down conversion. 

With the results of figure 5 in support, and utilising the 
time averaged energy expectation results of figure 4 as a 
guide, we then proceeded to compute the energy transfer, 
via the SQUID ring, from input field to an undriven out- 
put tank circuit at selected flux bias points close to and 
within an exchange region in 4. Our results are shown in 
figure 6 for the region of bias flux denoted using the time 
averaged energy expectation values at the top (which re- 
produces figure 4). As can be seen, outside the exchange 
region there is no evidence of energy conversion between 
the input microwave field and the output tank circuit. 
However, within the exchange region there is very sig- 
nificant conversion, reaching a maximum gain (w 70dB 
above background) at its centre. 

The results we have presented in this paper demon- 
strate that a input-output frequency down conversion of a 
factor of a couple of orders of magnitude can be explained 
through the simple model we have utilised. However, this 
does not constitute a theoretical upper limit, only a prac- 
tical constraint arising from the level of computational 
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FIG. 6: (top) reproduction of figure 4 illustrating the region of 
bias flux used to compute (bottom) the power gain in the tank 
circuit (with respect to $J tat = 0$o). Here, the microwave 
source frequency was 13GHz which converted into the output 
(tank circuit) oscillator mode of 130MHz. As can be seen, the 
peak power at the centre of the exchange region is ~ 70dB 
above the background from outside the exchange region. 



power available to us. We also note that by symmetry, 
and due to the fact that the SQUID ring couples differ- 
ent frequencies together, frequency up conversion via this 
same mechanism should also be possible. 



CONCLUSIONS 

In this paper we have demonstrated that the interac- 
tion of a quantum mechanical SQUID ring with classi- 
cal circuit environments is non-trivial. However, in the 
usual approach to the influence of a classical (and dis- 
sipative) circuit environment on the time evolution of 
a quantum mechanical SQUID ring [12], it is assumed 
that the environment can be modelled by a bath of lin- 
ear harmonic oscillators linearly coupled to the ring [29] . 
This need not be the case. The highly non-perturbative 
nature of the SQUID ring in the quantum regime (and 
other Josephson weak link based circuits) means that the 
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ring-environment interaction can be very non-linear and 
may lead to unexpected results. One of these is clearly 
the extreme frequency ratio conversion possible between 
classical fields via a quantum mechanical SQUID ring. 
Since the problem of environmental (dissipative, deco- 
hering) effects is so central to the successful implemen- 
tation of quantum technologies, it is our opinion that 
for non-linear devices such as SQUID rings the role of 
the environment has to be reconsidered within this non- 
perturbative (non-linear) context. 

On a purely experimental level the theoretical results 
generated in this paper indicate that quantum SQUID 
rings can be used for very large frequency ratio down 
conversion between classical fields. This seems to be sup- 
ported by experiment [17] and may prove to be of con- 
siderable practical significance, especially with the cur- 
rent interest in classical THz communications technolo- 
gies [30]. 
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